Optically engineering the topological properties of a spin Hall insulator 



Balazs Dora/'B Jerome Cayssol,^''^''* Ferene Simon/ and Roderieh Moessner^ 

^Department of Physics, Budapest University of Technology and Economics, Budafoki ut 8, 1111 Budapest, Hungary 
"^LOMA (UMR-5798), CNRS and University Bordeaux 1, F-33045 Talence, France 
^ Max-Planck-Institut fiir Physik komplexer Systeme, Nothnitzer Str. 38, 01187 Dresden, Germany 
'''Department of Physics, University of California, Berkeley, California 94720, USA 

(Dated: January 26, 2013) 

Time-periodic perturbations can be used to engineer topological properties of matter by altering 
the Floquet band structure. This is demonstrated for the helical edge state of a spin Hall insulator 
in the presence of monochromatic circularly polarized light. The inherent spin structure of the 
edge state is influenced by the Zeeman coupling and not by the orbital effect. The photocurrent 
(and the magnetization along the edge) develops a finite, helicity dependent expectation value and 
turns from dissipationless to dissipative with increasing radiation frequency, signalling a change in 
the topological properties. The connection with Thouless' charge pumping and non-equihbrium 
Zitterbewegung is discussed, together with possible experiments. 
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Introduction. Topologieal insulators (TIs) are a fo- 
cus of attention, not least due to their possible appli- 
cation in spintronics and quantum computation. They 
represent distinct states of matter with robust, topolog- 
ically protected conducting helical edge/surface states 
P, m. The importance of the spin-orbit interaction is 
reflected in their charge carriers having their spin locked 
to their momentum. In particular the two dimensional 
TI, namely the quantum spin Hall (QSH) state, has been 
predicted for a variety of systems including graphene Q , 
HgTc/CdTe ^ and InAs/GaSb Q quantum wells, lattice 
models and multi-component ultracold fermions in 
optical lattices 0-[lH. Nevertheless, in all of these, the 
gapless helical edge state originates from a subtle band 
inversion [l], 0] which requires careful Bloch band struc- 
ture en gine ering 0-111 as well as a high degree of sample 
control |12l4l4| . 

Bloch states and energy bands arise from spatially peri- 
odic Hamiltonians in condensed matter systems. Extend- 
ing the periodicity in the time domain through a time- 
periodic perturbation increases tunability of the Hamil- 
tonian: the temporal analogue of Bloch states (the Flo- 
quet states) can be manipulated via the periodicity and 
amplitude of the external perturbation. 

Recently, topological phases of periodically driven 
quantum Systems h^ve be;n characterized using Flo- 
quet theory, extending the time- independent topologi- 
cal classification |16l4l8| . Interestingly, novel topologi- 
cal edge states can be induced by shining electromag- 
netic radiation on a topologically trivial insulator, e.g. 
a non inverted HgTe/CdTe quantum well with no edge 
state in the static limit jl^. Besides, a time-dependent 
perturbation may also be harmful for the coherence of 
the edge/surface states of TIs by introducing dissipa- 
tion. It is therefore natural to investigate to what extent 
the steady state of a TI remains robust against time- 
dependent perturbations and how the electrical and mag- 




FIG. 1. (Color onhne) The quantum spin-Hall insulator (light 
yellow rectangle) with its helical edge state (counterpropagat- 
ing red/blue arrows) in a circularly polarized electromagnetic 
field with frequency lj and wave vector k. In the plane z — 
the rotating magnetic field B(t) — Bo{cosu}t,smijjt) is per- 
pendicular to the a" direction (vertical green arrows). A small 
tilting from the z axis does not influence our results. 



netic properties are altered. 

In this work, we consider the one-dimensional heli- 
cal edge state of a QSH insulator in a circularly polar- 
ized radiation field. When increasing the radiation fre- 
quency, the steady edge state is found to switch from 
a dissipationless charge pumping to a dissipative trans- 
port regime. We characterize those regimes by their dc 
and ac photocurrent responses and provide experimental 
proposals to measure them. Finally, we demonstrate that 
the photocurrent, the magnetization and the Zitterbewe- 
gung phenomenon are ruled by the very same unit vector, 
whose winding number determines a topological invari- 
ant for the system. Although, our predictions could be 
tested by experiments similar to those in graphene [20j 
and HgTe/CdTe quantum wells [2l|, they rely on a differ- 
ent coupling mechanism, that is Zeeman coupling rather 
than orbital coupling. 

Model.- We consider a QSH insulator located in the 
xy plane and radiated by a circularly polarized electro- 
magnetic field A{t) = AQ{cos{ujt — kz),sm{ujt — kz)) with 
wave- vector k and frequency w, whose sign determines 
the helicity of the polarization (Fig. [T]). The time- 



2 



dependent Hamiltonian of the QSH edge reads [22| 

Hit) ^ vpa'ip ^ eA,{t)) + g [a+e~"^' + h.c] , (1) 

where cr is the vector of Pauh matrices representing the 
physical spin of the electron, p the momentum along the 
one-dimensional channel, vp the Fermi velocity, and e 
the electron charge. The electric current operator|23 
is j = evpo'^. The circularly polarized radiation acts 
on both the orbital motion through the vector potential 

(t) = Aq cos Lut and on the electron spin through the 
Zeeman coupling g = ^ofiM-B^Oj 9eS being the effective g- 
factor and /xb the Bohr magneton. Nevertheless at high 
frequency, the orbital effect can be safely neglected ac- 
cording to a simple semi-classical argument (for a more 
rigorous treatment, see [l^). An electron travelling at 
the speed vp in an electric field Eq = Aqoj = cBq (c 
the speed of light) during a time l/w picks up an en- 
ergy vpeEo/uj from the vector potential which has to 
be compared to the smallest energy quantum it can ab- 
sorb, Huj (restoring original units). Hence in the regime 
VFeEo/uj <C hu, only the time-dependent Zeeman effect 
is effective, and in this respect, our effective Hamiltonian 
differs significantly from other studies on similar systems 
[i^ - lsil with dominant orbital effect. For typical param- 
eters [vp = 10^ m/s, laser power of 1 mW focused onto 
an area of 1 mm^, yielding Eq sa 600 V/m), this requires 
Lj ^ 0.5 THz, i.e. lasers operating in the far infrared or 
in the visible range. We also assume that huj is smaller 
than the bulk gap of the 2D insulator. 

Floquet states. - In order to study the steady state of the 
edge, we solve the time-dependent Schrodinger equation, 
idt^^pit) = H{t)'^p(t), with = in Eq. ©. Ap- 
plying Floquet theory 32, S^l, the solution of the time- 
dependent Schrodinger equation is written as 



1'p(0 = exp{-iEaip)t)<i>aip,t), 



(2) 



where Ea{p) is the Floquet quasienergy, and $a(p, i) = 
<^a{p,t + T) with T = 2tt /ui. From this, physically 
equivalent steady states can be created [32| by shift- 
ing the quasienergy En^aip) = Ea{p) -f nu and defin- 
ing '^n,a{p,t) ~ <f>Q(p, t) exp(ma;) where n is a relative 
integer. Then, the quasienergy and wavefunction are ob- 
tained as 



^c.{P,t) 



^J\ + a{vFP - w/2) 



2A V a cxpiiujt) X — a{vFP — w/2) 



(3) 



(4) 



where a = ±1, A = \/ g^ + [vfP — w/2)2. The quasiener- 
gies describe the opening of a gap of size g around to/ 2 
|23| . This photoinduced gap is located at momentum 
p = ui/2vF and stems from one-photon assisted processes. 
A given {p, t) describes the steady state where an ini- 
tial state with g = would evolve adiabatically if we 
switch on the magnetic field at t = — oo. 



We introduce the average energy [33f , which is used to 
identify the filled Floquet states [27[, in analogy to the 
stationary situation [23|, as 



^a(p) = *+(t)i?*p(<) = 



A + 



uj{vfP — '^/2) 
2A 



, (5) 



which is always single valued as opposed to the ladder of 
quasienergies En.aip). 

High and low frequency regimes.- It is natural to dis- 
tinguish high and low frequencies in terms of the ratio 
of the Zeeman coupling strength g and radiation fre- 
quency uj. More specifically the Floquet spectrum hap- 
pens to be gapped for \uj\ < Ag and gapless for |w| > Ag. 
In the low frequency regime, the bands are well sepa- 
rated by the photoinduced gap for any momentum, the 
(a = — l)-band being the fully occupied one. In con- 
trast, in the high frequency regime (|w| > Ag), the states 
of the (a = +l)-band become lower in energy than the 
ones of the (a = — l)-band within the momenta range 
cj- < vfP < w+ with 4a;± = cj ± ^ uP- — 16g'^ 23|. The 
band touching at \uj\ = Ag has a clear signature in the 
total energy which picks up a singular contribution as 



Etot = Es{g,uj) 



(kl - 4.g) 



3/2 



(6) 



for |w| > 4g, while E,{g,w) = po[Ljy4-gHn{2Wy/^/g)] 
is a smooth function of (uj — Ag). The lattice constant 
is denoted by a, po aj'KVF, and is a high energy 
cut-off. The exponent 3/2 appears also in the orbital 
contribution to the ground state energy of two dimen- 
sional Dirac fermions [s^]. 

Electromagnetic response and topological invariants. - 
The electromagnetic response of the QSH edge state is 
more easily detected than the singularity in the ground 
state energy Eq. (j6]). As a main signature, a dc pho- 
tocurrcnt (j) is generated along the edge whose direction 
is determined by the helicity of the circular polarization. 
Interestingly there is no accompanying ac-current in the 
absence of orbital coupling. Moreover the current opera- 
tor being j = ewFC^, such a dc current also corresponds 
to a steady state magnetization (cr^) along the edge. 

We have obtained the full dependence of the dc- 
photocurrent/steady state magnetization for any arbi- 
trary frequency within the bulk gap of the QSH insu- 
lator. Besides we demonstrate that the dc photocurrent 
is directly related to a topological property of the time- 
dependent Floquet state, that is the topological invari- 
ant: 

T 

C„ = i ^ / dtdo^At) ■ {dpdc.,p{t) X 3td„,p(t)) . (7) 

P 

This Chern number Ca, associated with the band a, is 
the winding number of the mapping, {p,t) — >■ da.p(t) = 



3 



i)<T$Q,(p, = a(gcosajt, (7sina;t,?;_Fp— a;/2)/A, be- 
tween the 1+1 dimensional extended Brillouin zone in 
{p, t) space and the unit sphere [l|, [l^ , the summation 
being taken over occupied bands. 

In the low frequency regime {\uj\ < Ag), the dc pho- 
tocurrent 



(j) 



(t) ^ ^, 
2tt '^p^ ' 27r' 



(8) 



is independent of the coupling strength the charge 
pumped within one cycle (T) being exactly the unit 
charge. This adiabatic pumped current has been con- 



sidered in Ref. [2^. As noticed by Thouless [35|, the 
integer charge pumped across a ID insulator in one pe- 
riod of an (adiabatic) cycle is a topological invariant that 
characterizes the cycle. Here this specific quantization of 
charge stems directly from the quantized Chern number 



Ca = - / dp 



as\g'n{ijj)v F 



2A3 



= -asign(a;), (9) 



the ground state being the filled a = — 1 band, yielding 
(j) = eC-/T. The dc current is therefore dissipationless, 
protected by a photoinduced gapj23j. 

At high frequency, \uj\ > Ag, the system undergoes a 
photoinduced band inversion and the Chern number 



Ca = -asign(w) 1 - 



(10) 



is no longer quantized (Fig. [2]) , reminiscent of the trans- 
fer of Chern number between equilibrium bands which 
touch. We note that the Chern number is continu- 
ous at the transition \uj\ ~ Ag and vanishes slowly as 
Ca = —a2g/uj for \u]\ ^ g. 

The corresponding dc photocurrcnt is 



ZTT 



s=±l 



(11) 



While the current still satisfies (j) ~ eC^/T for |a;| > Ag, 
it is dissipative and no longer quantized due to the band 
touching, in analogy with the photovoltaic Hall effect 
24| in graphene. The photocurrent approaches the fi- 
nite asymptotic value (j) = e5sign(w)/7r for ^ g, 
which can be regarded as the lowest order, linear response 
correction to the current in g, hence the weak-coupling 
regime (Fig. [5]). 

Proposal for a measurement setup.- In practise, the 
weak coupling regime g <^ |w| is usually realized. There 
a typical radiation field (magnetic field strength of the 
order of 10~^ — 10^^ T) yields a photocurrent of the 
order of 0.1 — 10 pA, depending on the effective g- 
factor values, which can be significantly enhanced (^cff ~ 
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FIG. 2. (Color online) The induced photocurrent (blue solid 
line) and the C_ Chern number (black dashed line) are shown 
as a function of the frequency lo. The photocurrent roughly 
behaves as (j) « esign(a;) min((;, |cj|/2)/7r and it is maximal 
at the transition \u}\ = Ag between the low and the high fre- 
quency regimes. The Chern number becomes non-quantized 
when band touching occurs a.t Ag = \uj\. 



20 — 50) for materials with strong spin-orbit coupling like 
HgTe/CdTe, InAs/GaSb, HgSe or BizSea. Such induced 
current can be detected in a contactless measurement. 
When the total area of the QSH insulator is exposed to 
the radiation field (i.e. the laser's spotsize is bigger than 
the area of the sample) , a circulating loop current flows 
around the sample as in Fig. [1] A perpendicular mag- 
netic field is induced according to the Biot-Savart law as 
Bind = /^o2\/2(j) /t^L with L the linear size of a square 
shaped sample (//q the vacuum permeability), staying 
roughly constant within the sample. For {j) = 1 pA and 
L = 1 micron, this gives Bind = 1 pT. This induced mag- 
netic field is within the detectability limit of an ac SQUID 
[2^. Finally standard 2 contacts measurement can also 
be used in order to detect the photocurrcnt. For a strip 
sample with laser's spot size bigger than the width but 
smaller than the length, backscattering is induced [ioj. 
which suppresses the photocurrcnt. 

So far we have considered the idealistic situation for the 
generation of the dc-photocurrent, namely zero chemical 
potential in the QSH edge modes, strictly vanishing or- 
bital effect and no inversion symmetry breaking. In the 
following we discuss how additional effects may influence 
the dc photocurrent (see also Q). 

Orbital effect and ac current response.- When the 
vector potential is taken into account (in the typical 
vpcAq/uj <C 1 regime), an ac current develops on top of 
the dc one as jj) ~ jdc.+ iac cos{ujt). We have solved Eq. 
(fT|) numerically |23| with the vector potential, and the 
results are shown in Fig. [3] The induced ac component 
stays always small compared to the dc one because the 
vector potential without the Zeeman term cannot cause 
spin-flip processes and is unable to generate any current. 
Indeed, the matrix element for optical transitions due to 
the vector potential is ^+{p,t)a^^-{p,t) = g/X. There- 
fore, the extended Kubo formula [27l| predicts the scaling 
of the ac component as jac ^ f Fe^n[23l|. 

Effect of the finite doping on the edge. So far we have 
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FIG. 3. (Color online) The induced dc (by the Zeeman term, 
solid line) ± the ac (by the vector potential, dashed lines) 
currents are plotted as a function of the vector potential, for 
several values of g, jac is always smaller than jdc for physically 
relevant parameters. The current behaves as (j) ~ jdc + 
jac cos{ujt). 



considered the optimal situation for photocurrent gener- 
ation, namely zero Fermi energy in the QSH state. In 
the case of a finite chemical potential, the dc photocur- 
rent vanishes gradually as we move away from half filling 
of the QSH edge states. The one-dimensional momen- 
tum acts as a polarizing effective magnetic field in Eq. 
([1]). For large momenta \p\ » (|a;|, \g\)/vF, this polariza- 
tion is so strong that the circularly polarized magnetic 
field hardly induces any magnetization, while close to 
the Dirac point (p ~ 0), the magnetic field represented 
by the momentum is very weak, and the circularly polar- 
ized field dominates over the momentum. The induced, 
helicity dependent magnetization originates from these 
states living close to p = u}/2vp, as indicated by the non- 
trivial Aharonov- Anandan phase in this region [23 1 . 



Inversion symmetry breaking and static magnetic field. 
We also consider the effect of a perturbation go<Jx in the 
Hamiltonian Eq. ([1} in order to mimic an eventual in- 
version symmetry breaking [22| and subsequent non- 
conservation (as in HgTe/CdTe quantum wells). This 
static Zeeman term opens a gap at p = whereas the dc- 
current is mainly built up from states near p = uj/2vp. 
Therefore the effect of inversion symmetry breaking on 
the dc-photocurrent is expected to be weak. Indeed we 
have checked that the dc current, and also the ac compo- 
nent (in presence of orbital effect), are almost identical 
to those of Fig. Elfor go < 5 [ii- 

Zitterbewegung. The trembling motion of the center 
of mass coordinate, is caused by interference between 
the positive and negative energy states (i.e. interband 
transitions) [H] . The topological invariant measures it 
indirectly through Aa,p{t) in Eq. ([7]). The position oper- 
ator satisfies dtx = vpcf^ =: v{t). Generalizing Ref. [s^ 



VF 



~ ([n o n + (1 — n o n) cos(2Ai) -I- sin(2At)nx] ctq) . 

(12) 



with n ^ da^plto), t ^ t — to, no n is the dyadic product 
and ctq is the spin configuration at < = io- 

Conclusion.- Radiation of a helical edge drives a tran- 
sition between nondissipative charge pumping at low fre- 
quency and a high frequency dissipative regime, reflected 
in the behaviour of the photocurrent. Note that for (neu- 
tral) atoms in optical traps, one can introduce a Zeeman 
term without any orbital counterpart, or fabricate chiral 
edge states with spin quantized parallel to the momen- 



tum 
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without any vector potential, the full transi- 
tion from dissipationless to dissipative charge pumping 
can then be followed. 
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SUPPLEMENTARY ONLINE MATERIAL FOR 
"OPTICALLY ENGINEERING THE 
TOPOLOGICAL PROPERTIES OF A SPIN HALL 
INSULATOR" 

DERIVATION OF THE ELECTRIC CURRENT 
OPERATOR 

The electric current operator can be obtained from the 
continuity equation, which, in one dimension, reads as 



dp{x) dja 



at 



dx 



= 0. 



(13) 



In our case, the charge density operator is p{x) = 
e{c\.{x)c^({x) + c^^{x)c]^{x)) , where Ccr{x) annihilates a par- 
ticle with spin a at position x. The time derivative can 
be evaluated from 



= i[H,p{x)] = -eu^(c|(x)ct(x) - c\{x)c^{x)) 



djx 
dx 



(14) 



From this, the current operator is obtained in second 
quantized form as jx = ew(c|c-|- — cjcj,), or in first quan- 
tized form as ~ eva^ . 



INFLUENCE OF THE VECTOR POTENTIAL 
AND OTHER TERMS 

Vector potential 

The effect of the vector potential can be inves- 
tigated more rigorously, after performing a unitary 
transformation as Vl'p(i) = Up{t)^p{t) with Up{t) = 
e'!cp(ia''vFeAQsin{ujt)/uj). As a result, the Hamiltonian 
changes (together with the —iUp{t)dtUp{t) term, coming 
from the time-dependent unitary transformation) to 

H ~ vpa^p + g [cr^ exp(— iwt — i2vi?eAo sm{u!t)/uj)+ 
+(T^ exp{iujt + i2vFeAosm{ujt)/uj)'^ . (15) 

The exponential term in the off-diagonal can be simplified 
using the Jacobi- Anger expansion as 

oo 

exp (i2zsin(a;t)) = Jm {2z) exp{imu!t), (16) 



z = vpeAo/uj, J„i{2z) is the mth Bessel function of 
the first kind. Higher harmonics in uit appear in tT+ 
and also cr^, triggering transitions between the eigen- 
states of tr^. This can be simplified in the limit of 
l^l ^ 1 using the expansion of the Bessel functions 
(J„i(2z) — > (sign(m)2:)l'"l/l'^l' fo^' small z), which trans- 
lates to vpeEo <g; huj"^ (upon reinserting original units). 
In this case, only the m = component needs to be 



taken into account, the higher harmonies can safely be 
neglected or treated perturbatively. 

In terms of the higher harmonics of the vector poten- 
tial, we can take them potential perturbatively into ac- 
count. The higher harmonics in ujt from Eq. ap- 
pearing after the unitary transformation, would dynam- 
ically open small gaps in the DOS, in addition to the 
main dynamical gaps from one-photon processes around 
±oj/2 with size g, at integer multiples of a;/2. Their size 
is estimated in the g lo limit as 



g f vpeEo 
\m\l \ up- 



|m 



g for TO 7^ 0, 



(17) 



corresponding to multi-photon excitations, similarly to 
graphene radiated with circularly polarized light (23. 27 1. 



While the Aq dynamical gap is generated solely by the 
circularly polarized magnetic fields, higher gaps are trig- 
gered by the vector potential term, though these are neg- 
ligible. 




FIG. 4. (Color online) The Floquet quasienergies are shown in 
the Floquet Brillouin zone (|-Bn,a(p)| < l^^^l/S) for g = O.SctJ. 
The black solid line/blue dots denote the analytical result 
from Eq. (3) in the main text/the numerical solution of Eq. 
(1) in the main text with A^; = 0, hardly distinguishable 
from each other. The red dots stem from the numerical solu- 
tion of Eq. (1) in the main text with z = vpeEo/ijj'^ — 0.2, 
the dynamical gaps are denoted by the green vertical arrows, 
showing good agreement with the prediction of Eq. (|17[1 as 
Ao = 0.3(^, Ai = O.OGoj, A2,3 are invisible on this scale. 
The vector potential is indeed neghgible for vpeEo <^ uJ^ , as 
conjectured using the Jacobi- Anger expansion, and only the 
circularly polarized magnetic field needs to be considered in 
this parameter range. 

To check the opening of additional dynamical gaps, 
we have evaluated the Floquet quasienergies of the 
edge state by the physically transparent, yet unortho- 
dox method of following the adiabatic time evolution 
of the initial eigenstates of Eq. (1) in the main text 
with Aq = g ^ 0, namely (0,1)'^ and (1,0)^. Both the 
vector potential and the Zeeman coupling are switched 
on at t — > — oo in an adiabatic manner, by attach- 
ing the exp(— i5|i|) factor {5 — >■ O"*") to them, and let- 
ting them evolve till t = 0. The resulting quasiener- 
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gies within the Floquet Brillouin zone [s^, HI] are deter- 
mined from the numerically evaluated wavefunction as 
E = u}\D.[^{t = 0)/*(t = T)]/2TTi, devised by Eq. (2) 
in the main text, shown in Fig. 21 As seen, the nu- 
merical procedure agrees convincingly with the analytic 
prediction of Eq. (3) in the main text. The main gap 
at vfP = ti^/S of size g from the Zeeman term domi- 
nates and the higher order dynamical gaps due to the 
vector potential are indeed negligible for vf^Eq/oj'^ <^ 1, 
in accordance with Eq. (|17p . Among the higher order 
dynamical gaps, Ai opens a tiny gap around the Dirac 
point E ~ 0. 

We have also solved our Floquet problem following the 
more conventional approach as in Ref. [stI by making use 
of the time periodicity of the Hamiltonian. There, by 
close analogy to Bloch states, the initially time depen- 
dent Schrodinger equation can be represented as a time 
independent matrix equation, whose various entries de- 
note the matrix elements between Floquet states. Then, 
the resulting eigenvalue problem is solved by truncating 
this matrix to a given size, corresponding to the allowed 
maximal photon excitations or the number of considered 
Floquet quasienergies. The obtained Floquet eigenener- 
gies and eigenfunctions are then used to calculate numer- 
ically the average energy to determine the proper filling, 
and then to evaluate the induced dc and ac component 
of the current. In general, the vector potential induces 
all higher harmonics of uj to the current as cos(na;i) with 
n integer. However, in the physically relevant case, when 
VpeAo ^ the cos(aj<) dominates over the other terms, 
and the resulting current reads as 

0') ~ jdc + jac COs(wt). (18) 

This is plotted in Fig. 3 of the main text, allowing for 80 
photons or equivalently 80 Floquet bands. We have also 
checked that the results do not depend on the number of 
considered Floquet bands in this range. Usually, jac ^ 
jdc in the (g^vpeAo) <ti uj regime, because the vector 
potential in itself cannot induce any current, only when 
the Zeeman term is included. 

In the regime of adiabatic charge p umping 
(u), VpeAo g), the Goldstone-Wilczck formula [22|, |38[ 
can be used to evaluate the current from Eq. as 

(j) = ^ (cj + 2vFeAo cos{ujt)) , (19) 

ZTT 

corroborating our numerical findings. Note that a finite 
g is essential to induce any current, although its explicit 
value drops out from the above expression. 

Inversion symmetry breaking 

We also consider the effect of an additional term in 
the Hamiltonian, which is a static Zeeman term in the x 
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FIG. 5. (Color online) The Floquet quasienergies are shown in 
the Floquet Brillouin zone in the presence of a static perpen- 
dicular Zeeman field, mimicking inversion symmetry break- 
ing. Top left: go = 0.2aj, g = Ao — 0; top right: go = 
VpeAo = 0.2aj, g = 0; bottom left: g = go = vpeAo = 0.2ci;; 
bottom right: g = go = 0.2a;. The main effect of the static 
Zeeman term is the opening of a gap ~ go around p = 0. 

direction as 

H' = goa\ (20) 

i.e. perpendicular to the spin quantization axis of the 
QSH edge state. A cr*' term would have identical effect. 
This can be thought of as mimicking inversion symme- 
try breaking^, since the edge state gap induced by the 
perpendicular static magnetic field is non-zero. Usually, 
the static in-plane field go is small and does not modify 
qualitatively the physics discussed so far. We have eval- 
uated numerically the Floquet bandstructure due to this 
additional static field in the presence and absence of cir- 
cularly polarized Zeeman and orbital terms. The main 
effect of go is to open a gap around p ~ of size go^ as 
can be checked in Fig. O 

In terms of the induced current, when g =■ 0, the ad- 
ditional static term can only induce a tiny ac current 
without any dc component, whose magnitude is compa- 
rable to that shown in Fig. 3 in the main text. When 
5o 9i which is the physically relevant regime, the in- 
duced current is almost identical to that in Fig. 3 in the 
main text, including both a dc and ac components. Only 
when go is comparable to the above picture is modi- 
fied and even in the absence of a vector potential, both 
dc and ac currents are induced by the interplay of the 
static and circularly polarized Zeeman fields. 

Two circularly polarized electromagnetic fields 

The case with two circularly polarized electromagnetic 
fields with different frequencies (wi and lij-2) can be con- 
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sidercd via the Hamiltonian 



H{t) = VF<7'-p + [a+ {gie 



— iujit 



526 



h.c. 



(21) 



neglecting the vector potential for simplicity. In the topo- 
logically protected region (a;i^2 ^ 51,2), when the peri- 
ods are commensurate to each other with lowest com- 
mon multiple T — ciTi = C2T2 (ci^2 integers), adiabatic 
charge pumping occurs. The dc photocurrent, which is 
the quantized charge adiabatically pumped through the 
system per cycle T.^s calculated from the Goldstone- 
Wilczek formula 13, 2^ as 



(j> = ^ ((^1 - W2)sign(gJ - gl) 



- LOl + CJ2) 



(22) 



giving (j) = ewi/27r eci/T for gi > 52 and (j) = 
ea;2/27r = ec2/T for 172 > ffi- For incommensurate 
frequencies, continuity suggests that this relation still 
holds [3^. In addition, the two frequencies also induce 
an ac current with lowest harmonics as 



3ac{t) = - 



65152(3? - 5l)(^^i - ^2) 



Mat + 9l) 



2\2 



COs[(wi — L02)t\. 



(23) 



For large frequencies (and small couplings), we can use 
linear response theory to obtain the average dc current 
as (j) = e [gisign(a;i) -I- g2sign(aj2)] /ti". 

This also explains quantitatively what happens to the 
Floquet edge state in a static magnetic field. By setting 
UJ2 = 0, 52 represents a static Zeeman term in the x 
direction as 5217^. From Eq. (P^ . the dc current remains 
unchanged for 52 < 5i with respect to its zero static field 
value as (j) = euijl-K^ and drops to zero for 52 > 5ii i-S- 
the static field destroys the adiabatic charge pumping 
when it becomes comparable to the circularly polarized 
component. 

In summary, the typical hierarchy of energy scales in 
a condensed matter realization of the QSH edge state 
is 5o < g < vpqAq < lo. In this case, as we have 
demonstrated, the induced dc current is mainly deter- 
mined by the circularly polarized Zeeman term, while an 
additional, smaller ac component cos{Lot) arises from 
the orbital effect. In cold atoms, the Zeeman term can 
be realized without the orbital counterpart, thus only a 
purely dc current in induced. 



DENSITY OF STATES 

The density of states (DOS, p{E)) is evaluated from 
the the overlap of two wavefunctions, one in which a par- 
ticle is created in a given state in the initial wavefunction 
and then it is evolved in time until and the other, where 
the initial wavefunction is evolved in time and then an 
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FIG. 6. (Color online) The quasienergies (left panel) and the 
corresponding density of states (right panel) are shown for 
g/uj = 0.2 (upper) and 0.8 (lower). The spectral weight in 
the spectrum increases from the blue and thin lines to red 
and thick ones. The black dashed hue denotes the average 
energy from Eq. (5) in the main text. For g = the DOS is 
constant. In the static case, a gap open of size g opens around 
zero energy. With increasing uj, the gaps are shifted to ±aj/2, 
carrying only half of the spectral weight. For g > uj/2, a clean 
gap remains around zero energy, but new gap edges show up 
at ±u]/2 ± g. For larger frequencies, the gap disappears, and 
a plateau develops around zero energy in the DOS at po/2. 



extra particle is added at t. The temporal Fourier trans- 
form of the time dependent overlap yields 

00 

g^{E) - i j dt cxp[iEt - tS]'^^{t)^p{0) - 


00 

/dt 
— cxp[i{E + Eait))t - tS]x 



X [A + a{vFP — w/2) + cxp{—iujt)(X — a{vFP — ^^72))] , 

(24) 

where i5 — )■ 0+. This gives (see Fig. 2 in the main text) 

(25) 



with po = a/'KVp^ a the lattice constant. Following the 



steps outlined in Refs. |2J, [251, and |27|, one obtains the 



very same expression for the DOS, shown in Fig. [6] to- 
gether with the Floquet and average energies. In addition 
to Ea{p), another branch with quasienergy Ea{p) — uj 
appears due to the spinor structure of the wavefunc- 
tion. The spectrum with positive or negative velocity 
{dpEa{p) > or dpEc,{p) < 0) hosts dominantly up or 
down spin electrons, respectively. Due to the circular 
magnetic field, these components are mixed and for small 
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FIG. 7. (Color online) The proposed experimental scheme for 
the detection of the photoinduced current. A linearly polar- 
ized radiation field is transformed to circularly polarized one 
by an electro-optical modulator (EOM)[43l, and the helicity 
of the polarization is changed periodically in time, with a fre- 
quency in the few 100 Hz range. The sample is denoted by the 
brick wall lattice filled plane. The direction of the photocur- 
rent (red and blue arrows within the plane of the sample) as 
well as the induced magnetic field (red and blue arrows per- 
pendicular to the plane of the sample) depend on the helicity 
of the polarization, Bind is measured with an ac SQUID. 

m, the up/down spin branch develops a weak down/up 
spin character, respectively. 

AHARONOV-ANANDAN PHASE 

Through the time dependent Zccman effect, the spin 
evolves cyclically with period T, and the wavefunction 
picks up the Aharonov-Anandan (39| phase: 



(26) 

For large momenta, this dynamical phase hardly changes 
with the parameters g and w, and is pinned to the trivial 
value 7 ~ — 27r8(— ap). By contrast, for p = w/2, it 
takes the non-trivial value — tt, signaling that states close 
to p ~ bjjlvp are the most influenced by the radiation 
field. 

POSSIBLE EXPERIMENTS 

Technically, a linearly polarized light can be run 
through an electro-optical modulator (EOM) , producing 



circularly polarized light (40|. as shown in Fig. [T) This 
can also change the helicity of the polarization periodi- 
cally in time, with a frequency in the few 100 Hz range 
(i.e. much smaller than that of the radiation field), result- 
ing in a sign change of the induced magnetization, i3i„d 
(pointing upwards vs. downwards with respect to the 
plane of the QSH insulator), facilitating the experimental 
observation. The induced magnetization stemming from 
the loop current can be detected by an ac SQUID (purple 
frame around the sample). By using the Biot-Savart law, 
the magnetic field profile along a cut through a square 
shaped sample is shown in Fig. [H which divides it into 
two identical rectangles. Qualitatively similar fields are 
induced along other directions. The divergence at the 
edges is cut-off by the finite spatial extension of the edge 
states. 
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FIG. 8. (Color online) The induced magnetic field is shown, 
calculated from the Biot-Savart law, staying rather flat within 
the sample (0 < a: < L). 



The induced current can also be detected directly using 
more conventional techniques. By irradiating only one 
edge of the sample and adding contacts to its ends, one 
could in principle measure the induced photocurrent, as 
shown in Fig. [S] 




FIG. 9. (Color online) Radiated single edge with contacts, 
designed the measure the photoinduced current directly. 



